Tripartite thermal correlations in an inhomogeneous spin-star system 
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We exploit the tripartite negativity to study the thermal correlations in a tripartite system, that 
is the three outer spins interacting with the central one in a spin-star system. We analyze the 
dependence of such correlations on the homogeneity of the interactions, starting from the case 
where central-outer spin interactions are identical and then focusing on the case where the three 
coupling constants are different. We single out some important differences between the negativity 
and the concurrence. 
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I. INTRODUCTION 



The notion of thermal entanglement relies on the 
amount of entanglement possessed by a physical system 
when it had undergone a thermalization process that has 
brought it in a thermal state In fact, in the presence 
of interaction between different parts of a compound sys- 
tem, even after a complete thermalization, the system 
can exhibit appreciable quantum correlations since the 
Hamiltonian cigenstatcs are, in general, entangled states. 
The establishment of a relation between temperature and 
entanglement has quickly brought to the idea of using the 
entanglement as an order parameter in quantum phase 
transitions Moreover, it has given a stronger stimu- 

lus of searching quantum correlations in the macroscopic 
world 



Thermal entanglement has been studied in connection 
with many possible applications and in different physi- 
cal systems: detailed analysis in spin chains described by 
Heisenberg models has been given Q as well as in atom- 
cavity systems Q and in simple molecular models Q. 
The quantumness of correlations has been singled out 
in thcrmalized spin chains (l0| also in connection with 
non local effects , and application of thermal entan- 
glement in optimal quantum teleportation protocols has 
been proposed [11]. 

Spin systems have been studied in depth not only in the 
neighbor-interaction configuration, leading to the Heisen- 
berg model, but also in star networks. This kind of sys- 
tems can be realized in many physical contexts, from 
Josephson Junctions [l| to trapped ions Q to solid 
state physics (l5j. In 2004, Hutton and Bose have ana- 



lyzed a physical system made of a central spin interacting 
with a set of N outer spins at zero temperature, bring- 
ing to light interesting properties which are immediately 
traceable back to the features of the ground state of this 
spin-star network fl6j . They have shown that the even- 
ness or oddness of the number of outer spins determines 
the law the entanglement amount scales with. The same 
spin configuration, in a simplified version involving three 



outer spins only, has been recently considered by Wan-Li 
et ai [ItJ in the special case where all the interactions 
between the central spin and the outer ones are identi- 
cal. Exploiting the high degree of symmetry of the sys- 
tem, they concentrate on appropriate concurrences [18[ 
to extract information about the presence of quantum 
correlations. However, the disclosure of quantum corre- 
lations in tripartite systems could require tools and con- 
cepts more adequate than the simple concurrences. Sabin 
and Garcia- Alcaine have contributed to solve this still 
very open problem fl9j introducing the notion of tripar- 
tite negativity [2(j, which is an effective tool to reveal the 
existence of quantum correlations traceable back to the 
impossibility of separating any of the three subsystems 
from the other two. 

In this paper we consider a spin-star system where 
three outer spins are coupled to the central one with 
different strengths, and analyze thermal entanglement 
in different situations. In the next section we present 
the model under scrutiny, which is characterized by 
anisotropic spin-spin interactions due to the absence of 
longitudinal z-z couplings. In the third section we con- 
sider the system in the homogeneous case, showing the 
presence of sharp changes in the amount of entangle- 
ment at zero temperature due to abrupt variations of 
the ground state. In the fourth section we analyze the 
inhomogeneous model, focusing on two types of inhomo- 
geneity. Finally, in the last section some comments and 
conclusive remarks are given. 



II. PHYSICAL SCENARIO 

System and Hamiltonian - In this section we present 
the spin star-system we have analyzed, which is pictured 
in Fig. [U Numbers 1,2,3 indicate the three outer spin 
1/2, while the central one is indicated with latin capital 
letter C. With ci,c 2 and c 3 we indicate the coupling 
constants of the central spin with the outer spin 1, 2 and 
3, respectively. Moreover, the whole system is immersed 
in a constant uniform magnetic field of modulus Bq and 
directed along the z-axis. The Hamiltonian model is then 
given by (with h = 1): 
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H = H!+H 2 , 
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FIG. 1: The physical system. Spins labeled 1,2 and 3 are 
coupled to the central (C) one. The magnetic field Bo is 
ortogonal to the plan where the four spins lie. 
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where a\ = |+>. <+| - a\ = |+) ? (-|, at = 

|— } i (+|, u>o is the unperturbed Bohr frequency, and Cj's 
are coupling constants. 

The term Hi, describes the interaction of the system 
with the magnetic field, while the second one, H2, arises 
from the dipole-like interactions between the central spin 
and the outer ones. The absence of z-z interaction re- 
flects a certain degree of anisotropy. 

Assuming that the system is in a thermal state, its 
density operator takes the form: 



P = 



Tr 



(4) 



so that p and H share the same eigenstates. The Hamil- 
tonian model given in Eq. ([1]) takes into account possi- 
ble inhomogeneities, nevertheless it is surely of interest 
studying it in the homogeneous case (ci = C2 = C3 = c) 
which turns out to be simpler from a mathematical point 
of view, paving the way to the study of the more compli- 
cated inhomogeneous model. Wan-Li et al. 17j have 
analyzed the homogeneous model taking advantage of 
the concurrence to quantify entanglement. In addition, 
they have assumed a more general anisotropic interac- 
tion. Since in real situations the construction of a sys- 
tem with perfectly homogeneous interactions could be 
quite difficult, in this paper we investigate the features 
of thermal quantum correlations in the presence of inho- 
mogeneity. 

Tripartite negativity - Remarking that the analysis 
is performed on the outer subsystem, which is made of 
three parts, we have the need to study quantum correla- 
tions in tripartite systems. There are various proposals of 
tripartite entanglement quantifiers (2l| - [2^ | and witnesses 
26|428j. the latter ones based on the key assertions in 
25| . Nevertheless, up to now and to the best of our 



knowledge, there is not a definitive answer to the request 
of quantifying tripartite entanglement For instance, 
it has been shown recently that three-tangle [29| could 
be an improper tool for such a purpose [30j. The situ- 
ation is much more complicated when we have to con- 
sider mixed states. Indeed, all the recipes mentioned in 
[21M24I [29j refer to pure states and fail when applied to 
non-pure states. Since we are studying thermal corre- 
lations, therefore dealing with highly mixed states, we 
need a tool to investigate quantum correlations in mixed 
states of tripartite systems. Again, to the best of our 
knowledge, the tripartite negativity introduced by Sabin 
and Garcia- Alcaine [2(| is a good quantity that allows to 
study quantum correlations in tripartite systems, even in 
the case of mixed states. Other quantities are very much 
related to the specific structure of the analyzed system, 
since observables that assume values in a certain range 
when the state is entangled are considered. Sabin and 
Garcia- Alcaine apply bipartite negativity to all the three 
possible bipartitions that can be singled out isolating one 
subsystem and considering the other two as a whole; then 
they consider the geometric mean of these three quanti- 
ties. According to the definition of negativity [HI, |32| . 
the partial negativity related to the bipartition I — (JK) 
(in which the two subsystems J and K are considered as 
a whole) is given by 



A/}. 



JK 



Ei 



(5) 



where ai(r TI ) is the z-th eigenvalue of t t/ , which is the 
partial transpose related to subsystem /, of the total 
(UK) density matrix. The parameter used to study cor- 
relations in tripartite systems is then 



M23 = v/M-23 A/2-13 A/3-1 



(6) 



which is the above mentioned tripartite negativity. In [20( 
the following properties have been proven: 

i) r separable or simply bi-separable => A/123 = 0; 

ii) Invariance of A/123 under LU operators; 

iii) Monotonicity of A/123 under LOCC operators; 

Though A/123 is not a sufficient condition to single out 
tripartite entanglement, it is an effective tool to study 
quantum correlations in tripartite systems. In fact, find- 
ing A/123 7^ implies that none of the three subsystems is 
separable, hence revealing correlations involving all the 
three subsystems. It is also important to note that neg- 
ativity Mi-jk does involve all the three spins, and es- 
tablishes the existence of a correlation between / and the 
entire couple made of J and K , which is a very different 
operation from tracing over one spin, say K, and eval- 
uating the degree of correlation between the other two, 
say / and J. 



3 



III. HOMOGENEOUS MODEL 

Let us consider our model in the special case Ci = 
c 2 = C3 = c that we have already addressed as the ho- 
mogeneous case. The complete diagonalization of this 
Hamiltonian model is given in Appendix A. In order to 
compute tripartite negativity we need first to find explicit 
form of the outer-spin density operator P123, tracing over 
the degrees of freedom of the central spin. 

It is quite clear that because of homogeneity, 
A/1-23 = A/2- 13 = A/3-12, so the geometric mean reduces 
to one of these three quantities. In Fig. [2] we plot 
A/123 as a function of temperature (T) and coupling 
constant (c). One can see some interesting features of 
A/123, f° r example its diminishing as T increases, and 
the presence of abrupt transitions at low temperature 
due to fast variations of the ground state. In Fig. 
El where A/123 (fes? 7 ) c = 6wo) is plotted, it is better 
shown this behavior which is due to the fact that the 
more temperature increases the more pem approaches 
a maximally mixed states (^I), making also P123 
maximally mixed. A more interesting trend of A/123 is 
observed at low temperature where abrupt variations 
are present, as shown in Fig. |4j At fc^T = O.OIcjo fast 
entanglement variations are well visible at c w 0.6o>o and 
c « 3.7wo. The reason for this occurrences is that in each 
of these points the system undergoes a level crossing 
involving its two lowest levels (ground and first excited 
states) . Let us call Ai , A 2 , A 3 the three plateaux of Fig. 
[2J corresponding to (c < 0.6wo) (0.6wo < c < 3.72uq) 
and (c > 3.72a;o), respectively. The entanglement 
in these areas is the same of the entanglement pos- 
sessed by the ground state. In A\ entanglement is 
evidently zero and, in fact, the ground state of the 
whole system is |g < -" 4l ' ) ) = |0000), corresponding to 

Pi23 l} = I 000 ) ( 000 l ( with thc notation \a c o x a 2 03)). In 
regions A 2 and A3 ground states of four spin system are 
\g( A ')) = (|0100) + |0010) + |0001) - V3 |1000)) /V6 
and \g {A3) ) ((|0011) + |0101) + |0110)) 

- (|1100) + |1010) + |1001)))/V6, and thc respective 
density matrix are p[ 23 } = (\W) (W\ + |000) (000Q /2 
and p[ A ^ = (\W) (W\ + \W)(W\j /2, where 

\W) = (| 100) + [010) + |001)) /VS and \W) = 
(|011) + |101) + |110))/\/3. 




FIG. 2: Tripartite negativity A/123 for the homogeneous model 
as a functions of coupling constant (c) and temperature 
(fcsT), both in units of ujq. Both the diminishing for increas- 
ing temperature and the low-temperature transitions with re- 
spect to c are well visible. 




FIG. 3: Tripartite negativity A123 of the homogeneous model 
versus fcsT, at c = 10u;o. Temperature is in units of cjq. In 
the area where ksT < ujq/2 we observe that the negativity is 
almost constant. 



IV. INHOMOGENEOUS MODEL 

Once we have analyzed thc homogeneous model, we 
turn to the study of the lack of homogeneity in the three 
interactions. For the sake of simplicity we shall concen- 
trate on two kinds of inhomogeneity that we call type A 
(ci = c, c 2 = ci, C3 = c) and type B (ci = c, c 2 = ex, 
C3 = cx 2 ), in both of which the inhomogeneity is char- 
acterized by a suitable dimcnsionless real parameter x. 
In both cases we study the thermal quantum correla- 



FIG. 4: Tripartite negativity A123 of the homogeneous model 
versus coupling constants (in units of ojo), at fcsT = 0.01o;o. 
Abrupt variations at c ~ 0.6wo and c « 3.7ojo are well 
visible. Three plateaux, A\ (c < Q.Qloq),A2 (0.6ojo < c < 
3.72a;o), A3 (c > 3.72tJo) are quite evident. 



tions of the tripartite outer subsystem as a function of 
temperature (T), coupling parameter (c), and parame- 
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tcr of inhomogcneity (x). It is worth noting that in the 
homogeneous model the structure of eigenstates of H is 
independent on the coupling constant (c) , so that the fea- 
tures of the correlations are essentially determined by the 
structure of the eigenvalues only. Instead, in the inho- 
mogeneous case the eigenvalues are c— and x— dependent 
and the eigenstates depend on x, therefore the thermal 
correlations are affected both by the crossing of levels and 
by the modification of the structure of the eigenstates. 
We will see that while the change of the eigenstates pro- 
duces smooth changes of A/123, on the other hand a level 
crossing, especially at low temperature, can produce a 
very sharp variation of the negativity. 



A. Model with inhomogeneity of type A 

In spite of the lack of some symmetry, the inhomo- 
geneous model of type A is still easily solvable and its 
diagonalization is given in Appendix A. 




FIG. 5: Negativity A/123 versus temperature (fcsT, in units 
of uq) and inhomogeneity parameter (a;), at c = 6ojo, for 
the inhomogeneous model of type A. The trend with respect 
to temperature is decreasing. The dependence of A/123 on 
the inhomogeneity parameter, at low temperature, exhibits 
different trends: initially, smooth changes, then, fast falling 
around x ~ 5.5 and zero value from there on. 



Since one can see that in the Area A3 of Fig. @] neg- 
ativity assumes its greatest value, we start the analy- 
sis of the inhomogeneous model of type A showing in 
Fig. [5] the complete dependence of A/123 on temperature 
and inhomogeneity parameter for c = 6wo (this specific 
value of c belongs indeed to the area A3 of the homo- 
geneous model). It is quite obvious that the physical 
reason of the decreasing of A/123 with respect to T is the 
same as for the homogeneous case, i.e. the high degree 
of mixedness. At low temperature, in the < x < 1 
area, the smaller is x, the closer to zero is the negativ- 
ity. In fact, when x is quite smaller than 1, C2 = c x 
is much smaller than c\ — C3 = c, and hence, roughly 
speaking, spin 2 can be considered almost decoupled, 



[x + V8 + x 2 ) and 
+ x 2 + l)u>o, that causes a 



and then separable. Moreover, at xq « 0.43 this neg- 
ligibility of the coupling between spin C and spin 2 is 
concomitant to a level crossing between the two lowest 
energy levels, E^{c = 6u>o, x) = — 3wo (x + y/8 
E^(c = 6ujo,x) = — (6 \/2 + x 
sharp change at x = xq , where E^ = E± . 

Still, at low temperature, for x quite greater than 1, 
A/123 decreases as x becomes increasingly grater than 
1, because the coupling with spin 2 becomes stronger 
than the other two, making spin 1 and 3 almost sepa- 
rable from the spin 2. Around x ~ 5.5 a level crossing 
occurs. Because of the diminishing of A/123 for x ^> 1 
and x <C 1, there must be a maximum in the interme- 
diate region. Uncxpcctfully, this maximum is reached at 
x = xm ~ 2.46, instead of x = 1. We think this is an 
important result because, even though there is no big 
difference in the values of negativity in a wide contour 
of the maximum (let's say from x = 0.5 to x — 5.5). 
In fact, it is conceptually important that the maximum 
of correlations is not reached in the homogeneous case, 
where the high symmetry of the system could lead to the 
idea of stronger correlations between its parts. 

Evaluating the eigenvalues of Hamiltonian one 
can find that the ground state in the nearby of 
this maximum 



is 1^2); which becomes 



0.2073(|0011) 
0.6435(|1010) 
values of A/123 



*M / 

0110)) - 



|1100)) + 0.2073(|1001) 
1 0101)) for x = xm- In this region, the 
and in particular the appearance of a 
maximum arc essentially determined by the dependence 
of the structure of the state | 'F^ ) on the inhomogeneity 
parameter. The relevant density operator is a mixture 
of two Werner-like states, and can be written in the 
following way: 



P 

with 



■(x) 



- \wi{x)) { Wl {x)\ + - \w 2 (x)) (w 2 (x)\ , (7) 



K(x)) 

\w 2 (x)) 
where 



H 1011) 



|110) 



V8 



H 1 100) + |001) 



V8 



|101) (8) 



|010) (9) 



N = (2+ ((V8 + x 2 - x)/2j 



-1/2 



(10) 



The relevant negativity can be evaluated and it turns 
out to be: 



AV"(z)) = 



1 



1 



2 5 /3 ^ (g + ^2)3/2 



+ V8 + x 2 - yj2(2Q + x{x + V8 + x 2 )) 



-3\/8 + x 2 + 8v / 8 + x 2 
(|Ai(x)| + |A 2 ( a; )| + |A3( a; )|)] 2 } 



1/3 



(11) 



1 2 3 4 5 6 7 



FIG. 6: Comparison of the Negativity Af(p 2 ~(x)) (solid blue 
line) with the Negativity of the thermal state (dashed red 
line) evaluated for c = 6a>o and fcsT = O.OIojo, both vs x. In 
the region 0.5 < x < 5.5 (where / i s the ground state) the 
two curves coincide. 



with Aj(x)'s the roots of the algebraic equation in A: 

(-10a; + x 3 + (2 + a; 2 V8 + a; 2 ) 
+ (-32a; - Ax 3 - 48^8 + x 2 - Ax 2 y / 8 + a; 2 ) A 
+ (l28x + 16a; 3 - 384^8 + x 2 - 48a; V8 + a; 2 ) A 2 
+ (l024^8 + x 2 + 128a; 2 ^8 + a; 2 ) A 3 = 0, (12) 

where the coefficients, and hence the solutions, depend on 
x. The behavior of Af(p 2 ~ (x)) is shown in fig. [6j where 
the negativity of the thermal state for c = 6ujq and TsO 
versus x has also been plotted. It is well visible that the 
behavior of N{p 2 ~(x)) perfectly reproduces the edge of 
the top visible in fig. [5]in the low temperature region and 
for x belonging to the region where j^^) is the ground 
state of the system. 

It could be interesting to compare the behavior of the 
tripartite negativity with the values of the three concur- 
rences related to the three couples of spins that can be 
extracted tracing over one of them: C23 = C(Ti*ip 2 ~), 
C12 = C(Ti 3 p 2 -), C13 = C(Tr 2/ 9 2 -). Figure [7] shows that 
around x = 0.5 and x = 5.5 there are abrupt changes in 
the values of the concurrences (in the same points where 
the tripartite negativity exhibits the same behavior), due 
to rapid changes of the ground state. It is worth noting 
that in the region 3.5 — 5.5 the concurrences C12 and C23 
are vanishing while the tripartite negativity (and hence 
all the three relevant bipartite Negativity functions) is 
not. This apparent contradiction reflects the very dif- 
ferent meanings of these quantities. Indeed, for example, 
while the concurrence C12 measures the entanglement be- 
tween spins 1 and 2, without considering the spin 3 that 
has been preliminarily traced over, the negativity A/1-23 
takes into account the correlations between spin 1 and the 
whole couple made of the two spins 2 and 3. Therefore, 
in that region, we can talk about correlations between 
spin 1 and the couple (2,3) without correlations between 
1 and any of the other two spins, after the third has been 



FIG. 7: The three concurrences C12 = C23 (dashed blue line) 
and C13 (solid red line) vs x. Abrupt changes are very well 
visible around x = 0.5 and x — 5.5, where the negativity 
exhibits the same feature. For x = 1 the three concurrences 
assume the same value (« 0.33). 




FIG. 8: Negativity A/123 as a functions of temperature (fcsT) 
and coupling parameter (c) at x = 3, for the inhomogeneous 
model of type A, both ksT and c are in units of ujq. With 
respect to c, at low temperature, there are abrupt changes 
traceable back to the occurrence of level crossings. 



traced over. We think this is an important result, since 
it underlines the difference between Af and C . 

Figure [5] shows the complete dependence of A/123 on 
temperature and coupling parameter, for x = 3. We ob- 
served (by performing other plots, which we are not re- 
porting here) that the qualitative behavior of A/123 is in- 
dependent on the specific value of x, provided it is neither 
too small nor too big. It is quite evident that the trend 
with respect to temperature is the usual one. At low 
temperature the dependence of A/123 on c exhibits only 
the step trend due to level crossings. The highest step 
involves the energy levels E^(c,x = 3) = — cvTT — cjq 
and E^ic.x = 3) = — (3 + VT7) c/2 (see Appendix A), 
and it is quite easy to find that the crossing occurs at 
c t w 4.08wo where E± = while E± < E^ for c < c t 
and E2 < E^ for c > c*. 

To conclude the analysis of the inhomogeneous model 
of type A, in Fig. |H]we have plotted the complete depen- 
dence of A/123 on coupling and inhomogeneity parame- 
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are all present also in this model of inhomogencity. Since 
for type B we do not have an explicit analytical diago- 
nalization of the Hamiltonian, we also miss an analytical 
expression of the state of the system in the nearby of 
the maximum. In fact, we have only numerical results, 
according to which, for example, for T = O.Olwo/fcB 
and c = 6 uiq the maximum of negativity occurs at 
x = xm w 1.91, where the ground state of the system 



is approximately: 



0.292(|1100) - |0011)) 



0.471 (|1001) - |0110)) + 0.439 (|1010) - |0101)). 



FIG. 9: Negativity A/i 23, at low temperature (ksT = O.OIcjo), 
versus coupling parameter (c, in units of too) and inhomogene- 
ity parameter (x), for the inhomogeneous model of type A. 
Both smooth decreasing and fast variations occur with respect 
to x. Dependence on c exhibits a stepwise behavior. 



tcrs, fixing the value of temperature at fc^T = O.Olwo- 
The dependence on the coupling parameter is qualita- 
tively the same as for the previous case (see Fig. 
Indeed, we observe a step trend caused by a level cross- 
ing between the two lowest energy levels. As a function 
of x, A/123 exhibits both smooth decreasing, due to in- 
creasing separability of spins 1 and 3, and fast variations 
traceable back to a level crossing amongst the two lowest 
energy levels and . In fact, we verified that the 
big step observed in Fig. [3] corresponds to the locus of 
points of the plane (c, x) in which the two lowest levels 
have the same energy. Thanks to the simple form of the 
eigenvalues of our model, one can find a simple analytical 
expression of this curve: 



c(x) 



x 2 +4 + a(x) 1 x 2 + 5 + a(x) 



Ax 



where a{x) is 



a(x) = VT6 + 10x 2 + x 4 



(13) 



(14) 



B. Model with inhomogeneity of type B 

The second type of inhomogeneity we decided to study 
is characterized by the coupling constants c\ = c, 02 = cx 
and C3 = c x 2 . Also for this model we have considered 
the dependence on T, c and x. Surprisingly, what we 
have found is that the behavior of negativity in the pres- 
ence of inhomogeneity of type B is not different from 
what we have found in the presence of inhomogeneity of 
type A. The diminishing for increasing temperature, the 
abrupt changes at low temperature in connection with 
level crossings and the presence of a maximum for x ^= 1 



V. CONCLUSIONS 

In this paper we have analyzed the thermal correla- 
tions in a spin-star system consisting of a central spin 
interacting with three outer ones, all immersed in a mag- 
netic field. In the absence of a tripartite entanglement 
measure, we have decided to exploit the tripartite nega- 
tivity, which is at least able to put in evidence the lack of 
separability and simple bi-separability of a mixed quan- 
tum state of a tripartite system. 

We started from considering the homogeneous model in 
which the three peripheral spins interact in the same way 
with the central one, showing the appearance of a signifi- 
cant degree of inseparability revealed by a non vanishing 
tripartite negativity A/123. Typical behavior of thermal 
correlations reflects the behavior of A/123, which decreases 
to zero-value as temperature increases and, instead, ex- 
hibits both appreciable values and abrupt changes at very 
low temperature. These occurrences show in a very clear 
way the role of thermal entanglement mediator played by 
the central spin. Further we have considered the pres- 
ence of inhomogeneity, meant as differences in the cou- 
pling constants between the outer spins and the central 
one, focusing on two special cases. In the first case (in- 
homogeneity of type A) two spins interact in the same 
way with the central one and the third has a different 
coupling constant. In the second case (inhomogeneity of 
type B) all the three outer spins are characterized by dif- 
ferent coupling strengths. Though some differences are 
visible, the qualitative behaviors of A/123 are quite simi- 
lar for the two inhomogeneous models and even for the 
homogeneous one. This suggests the idea that the ther- 
mal quantum correlations mediated by the central spin 
are not very much damaged by a certain lack of homo- 
geneity that could characterize a more realistic situation, 
provided the degree of inhomogeneity is not high. 

A remarkable point is that we have singled out the 
presence of maxima of tripartite negativity correspond- 
ing to inhomogeneous models, i.e. for the values of the 
inhomogeneous parameter quite larger than 1. Though 
the differences between a maximum value of negativity 
and its values in the nearby (even up to the homogeneous 
case, i.e. up to x = 1) is not very big, we think that it 
is conceptually important the fact that a higher degree 
of symmetry in the system does not guarantee a higher 
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degree of correlations between all its parts. 

Another important point is that we have found some 
regions of the parameter space where it happens that 
two concurrences, say C\i and C23, are vanishing, mean- 
ing that there is no entanglement between 1 and 2 and 
between 2 and 3, while the relevant negativity A/2-13 is 
non-vanishing, meaning that there is a correlation be- 
tween 2 and the couple made of 1 and 3. This fact points 
out in a very transparent way the different meanings of 
concurrence and negativity. 



Appendix A 



In this appendix we give eigenvalues and eigenvectors 
of the Hamiltonian of the inhomogeneous model of type 
A (ci = c, C2 = c x, C3 = c), as functions of x > 0, c > 0, 
ujq. The special case x = 1 gives the solutions for the 
homogeneous model. 



Energies 



±c x , 



= ± 



E7 = CJQ , 



Eg 
Eg 



-2w 
2w , 



\ + x 2 )"- 



c(2 + x 2 



-w 



- w 



Eigcnstatcs 

= I [(10011) ± 11100)) - (|0110) ± 11001))] , 

1*3 ) = -k Kl 0011 > ± I 1100 )) + (l 011 °) ± I 1001 )) 

+ ^+f~* (|0101) ±|1010))1 , 

|*3 ) = XT [(l° 011 > ± l 1100 » + (l° 110 > ± l 1001 » 



(|0101) ± |1010>) 



1*4) = & [V2 



a; 2 |0111) 
±(|1011) + a; 1 1101) + 1 1110))] 
[(|0100) +x |0010) - 



1 

K 2 



\0 

±^f{2 + x 2 ) |1000) 
*e ) = jfe [s loooi) + 10010) 
*6 B > = Tifc, (10010) 

«tf> = £[*|ioii> + 
* 8 ) = loooo) , 

*9> = 11111) , 



-x|0001)) , 

liioi) - (i 

-x |1011)) , 



|0001)) 

(h+x) 10100)] , 

) lino)] , 



r 



where K\, K2 and K3 are: 



A1 2 =4 + 2 



V8 + x 2 ~ x 



A3 2 = -s + 3 + . 

or 



Ao 2 



2 (2 + x 2 ) , 
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